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ABSTRACT 
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We have imaged two normal, non-coronal, infrared-bright K giants, a Tau and a Boo, 
in the 1.4-mm and 2.8-mm continuum using BIMA. These stars have been used as im- 
portant absolute calibrators for several infrared (IR) satellites. Our goals are: (1) to 
establish whether these stars radiate as simple photospheres or possess long-wavelength 
chromospheres; and (2) to make a connection between millimeter-wave and far-infrared 
(FIR) absolute flux calibrations. To accomplish these goals we also present ISO Long 
Wavelength Spectrometer (LWS) measurements of both these K giants. The FIR and 
millimeter continuum radiation is produced in the vicinity of the temperature mini- 
mum in a Tau and a Boo. We find that current photospheric models predict fluxes 
in reasonable agreement with those observed for wavelengths which sample the upper 
photosphere, namely <125 /xm in a Tau and a Boo. We clearly detect chromospheric 
radiation from both stars by 2.8 mm (by 1.4 mm in the case of a Boo). Only additional 
observations can determine precisely where beyond 125 /xm the purely radiative models 
fail. Until then, purely radiative models for these stars can only be used with confidence 
for calibration purposes below 125 //m. 

Subject headings: infrared: radiation — infrared: stars — stars: atmospheres — stars: 
chromospheres — stars: individual (a Tau, a Boo) 

1. Introduction 

We had two distinct objectives in undertaking this work, dominantly to investigate the merit 
and reliability of K giant stars as absolute FIR and millimeter calibration sources, but secondarily 
to determine just how shallow in depth simple radiative equilibrium cool star model atmospheres 
might be trusted. 

We have chosen a Tau and a Boo for this initial study because these stars, being both bright 
and nearby, are easily the two most completely and thoroughly studied red giants. This makes them 
the ideal test cases for our purposes. It should be noted that they both lie on the cool, non-coronal 
side of the "Linsky-Haisch dividing line" (Linsky k, Haisch 1979) which separates coronal from 
non-coronal giants, although Ayres et al. (1997) have found evidence for buried high temperature 
emission in the case of a Boo. 

1.1. Infrared calibration stars 

Cohen and colleagues have presented a self-consistent context for absolute calibration in the 
IR, ideal for use by ground-based, airborne, and spaceborne spectrometers and radiometers (Cohen 
et al. 1992a,b,c, 1995; Cohen & Davies 1995; Cohen et al. 1996a,b, 1999). Their approach is 
based upon a pair of absolutely calibrated, IR-customized models of Vega and Sirius calculated 
by Kurucz. These efforts have furnished the absolute stellar spectrum of Sirius that underpins 
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COBE/DIRBE bands 1 - 5 (Mitchell et al. 1996; Cohen 1998), the spectra of KO-MOIIIs for 

the on-orbit calibration of the Near- and Mid-Infrared Spectrometers on the joint NASA/ISAS 
Infrared Telescope in Space (Murakami et al. 1997), likewise for ESA's Infrared Space Observatory 
(ISO) (Kessler et al. 1996) instruments (e.g., the Short Wavelength Spectrometer) (Schaeidt et 
al. 1996), for the Spatial InfraRed Imaging Telescope (SPIRIT-III) on the US Midcourse Space 
Experiment (MSX) (Mih et al. 1994), and for 2MASS (Cohen et al. 2003b). An extension of this 
work to substantially fainter stars has provided calibrators for the Spitzer Space Telescope (Cohen 
et al. 2003a). Most recently, Price et al. (2004) have absolutely validated the two calibrated model 
spectra, and the empirical spectra of the brighter K giants such as a Tau and a Boo, in six bands 
in the mid- infrared (MIR), using precision measurements from the MSX. 

The absolute calibration of ISOPHOT (Lemke et al. 1996; Schulz et al. 2002) was supported 
by planets, asteroids, and stars, ranked in order of decreasing flux density (Muller & Lagerros 
1998; Schulz et al. 2002). The stellar calibrators (Cohen et al. 1996b) necessitated long- wavelength 
extrapolations, as far as 300 fim, of the observed (1.2-35 jum) absolute stellar spectra by model 
atmosphere spectra (computed by one of us [DEC]: Cohen et al. (1996b)). We note that the bright, 
secondary reference stars that have been selected for IR calibration purposes are precisely those 
"quiet" giants in which it was suggested that radiatively-cooled regions very largely dominate the 
stellar surfaces so that single component atmospheric models were most likely to be valid. Examples 
of these are a Boo, a Hya, a Tau, and 7 Dra (Wiedemann et al. 1994). Consequently, one of our 
objectives is to establish whether these stars radiate as predicted, thereby validating the scheme 
used for ISO at the longest wavelengths, or possess long- wavelength chromospheres. 



1.2. Linking infrared and millimetric calibration 

The ISO Long Wavelength Spectrometer (LWS) (Clegg et al. 1996; Swinyard et al. 1996) used 
Uranus as a calibration standard but also observed Neptune, and Mars when it became available. 
Although Uranus has remained the LWS primary calibrator, the Mars spectrum has been shown 
to be fully consistent with this Uranus framework (Sidher et al. 2003). In the case of Neptune, a 
new model based on the LWS data has been adopted (Orton et al. 2000) . LWS was able to observe 
some stellar calibrators (Lim et al. 1997), although these observations were difficult to make due 
to poor signal-to-noise in the relatively narrow LWS bandwidths. However, the stellar spectra did 
demonstrate consistency with the Uranus-Mars framework, in the LWS 45-170 fim region. In this 
paper we present the totality of results from the ISO LWS spectra obtained on a Tau and a Boo. 

The long-wavclcngth stellar spectral extrapolations arc based upon the assumption that these 
calibrator stars arc not attended by extensive chromospheres. Radiometric testing with ISOPHOT 
(Muller & Lagerros 1998; Schulz et al. 2002) suggests that these extrapolations were meaningful, 
i.e. consistent with calibrations based on asteroids and planets. It would be of great interest, both 
scientific and pragmatic, to determine independently whether these same stellar calibrators radiate 
as predicted in the millimeter (mm) domain. The basis for mm-wave absolute flux calibration is 
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Mars. Therefore, we observed the stellar calibrators against Mars, or against secondary microwave 
standards that were ultimately traceable to a comparison with Mars, Uranus, or Neptune. By this 
means we hoped to provide a link between the reference flux density scales of planets and stars and 
to unify flux calibration across the IR-mm regime. Early efforts in this direction at Berkeley were 
described by Wrixon et al. (1971). Gibson et al. (2005) and Welch & Gibson (2005) have continued 
this radio absolute calibration work. 

1.3. Cool stellar atmospheres 

Our secondary goal is to probe the structure of the outer atmospheres in K and, eventually, 
M giants. Vernazza et al. (1976) demonstrated in a classic paper how observations of the solar 
FIR continua, atomic lines, and UV continua could be brought together to construct a single- 
component, time-independent mean solar temperature structure. FIR continuum fluxes were of 
crucial importance in helping to establish this mean structure from the upper photosphere, through 
the temperature minimum region, and into the chromosphere. Using the FIR continuum to probe 
the shallow depths in solar and cooler stars is possible because the primary IR continuous opacity, 
H~ free- free, increases as A^, pushing the depth of continuum formation into the outermost stellar 
layers. Until now these layers in cool stars have been accessible only through non-LTE analyses of 
the cores of strong atomic and molecular lines (e.g., Ayres &; Linsky (1975), Kelch et al. (1978), 
McMurry (1999), Wiedemann ct al. (1994)). 

While cliromospheric models have become more complex in the years since the Vernazza et al. 
(1976) study, the FIR and millimeter fluxes still have the potential to be valuable diagnostics of 
the shallow stellar layers. In accurately determining the FIR and millimeter fluxes of selected red 
giants, we hope to provide useful information for testing models of the cool star stellar photospheres 
and chromospheres. We are particularly interested in determining the location of the transition 
from observed fluxes, reproducible by purely radiative equilibrium models, to fluxes that require 
invoking a non-radiatively heated chromosphere. In addition, we wish to trace the observed run 
of decreasing fluxes with increasing wavelength from the upper photosphere and determine how 
closely current stellar model atmosphere calculations can match this run. This information on the 
location of the temperature reversal and the quality of fit to the fluxes originating in the upper 
photosphere are vital in defining the wavelength range for which simple radiative models may be 
used for FIR calibrations. 

2. The observations 

Table 1 summarizes our various sets of mm observations, listing array configurations and fre- 
quencies used, integration times on the stars, phase and flux calibrators (chosen to be as close to 
each star's direction in the sky as possible), and the primary or secondary mm-calibrator observed 
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with each set of stellar observations. A few measurements were taken in the B-array, with which 
planetary calibrators like Mars were partially resolved, necessitating an extra step in the calibration 
process to calculate the flux expected on the shortest baseline pairs (see below). Consequently, we 
preferred measurements in C- or D-array configurations. Angular diameters of these stars in the 
optical to MIR regions are of order 20 mas, rendering them unresolved targets in all these config- 
urations, as were the local calibrators chosen. On the basis of estimates of the BIMA performance 
at 1 and 3 mm, our observing scripts called for typical tracks of several hours' duration, of which 
between roughly two and four hours were spent on the stars themselves. 

This was a continuum experiment, and we selected observing frequencies from considerations 
of the expected stellar spectra, receiver performance, and atmospheric transmission. Whenever 
possible we included measurements of Mars with each set of stellar and local calibration source 
observations, in an attempt to unify the calibration. For those few data sets taken in an array 
configuration that resolved Mars, we calculated visibilities as a function of baseline (in kiloA) and 
applied these partial-resolution factors to the Mars model (Wright 1976) flux densities using our 
shortest baseline pairs. 

Sometimes no planet was visible close in time to our observing tracks. Then we used MWC 349 
when it was available. Hat Creek has a long series of observations of this star, it is relatively bright, 
has shown no statistically significant variability in our archives, and has been modeled in detail by 
Dreher & Welch (1983) so that flux extrapolation and interpolation to our frequencies is viable. 
Further, Welch &; Gibson (2005) have obtained its mm flux densities to a precision of 1% absolute. 
For our final observations, in December 2000, we used W30H as reference for our a Tau track, 
again because of its relative proximity to this star, long-term temporal stability of signal, and 
the available multifrequency Hat Creek archives which enable interpolation between known flux 
densities at other frequencies to our own. Wc treated W30H as having a canonical Hll region 
spectrum and included the frequency dependence of the Gaunt factor. Other secondary references 
such as 1415+133 and 0530+135 were likewise observed when no other obvious flux reference 
presented itself. Independent observations of these were extracted from the Hat Creek archives, 
preferably with respect to a primary or acceptable secondary standard. In such circumstances. 
Table 1 carries the primary standard to which the flux pedigree of such a secondary is traceable. 

Upper and lower sidebands were processed independently, maps made of the stars at each 
frequency, and peak flux densities extracted (see the final column of Table 1). Sometimes, weather 
conditions were adequate in the lower sideband but led to unacceptably large noise in the upper 
sideband. In such cases we eliminated the noisy datascts which were readily recognized by the 
presence of artifacts in the maps, such as global striping, or large negative-going blobs. To estimate 
the uncertainties in the tabulated peak flux densities derived from any map, we used the standard 
deviation calculated from substantial regions of the map that excluded the target star itself. To 
combine the simultaneous upper and lower sideband maps at one epoch, we used inverse-variance 
weighting of the two maps based upon these noise estimators, relying upon the MIRIAD task MATHS. 
Likewise, we combined results across epochs. 
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3. Self-calibration of the 1-mm data 

Phase calibration of interferometer data relies on transferring phases accurately from a point 
source (usually a QSO) to the target source. The accuracy of the calibration depends on the 
distance between source and calibrator, baseline errors, atmospheric phase stability, and the signal- 
to-noise ratio (SNR) in the phase measurement of the calibrator. All soTirccs of phase error result 
in decorrclation of the calibrated image; i.e., a lower amplitude in the final map compared with 
perfectly calibrated data. Self-calibration offers a way to eliminate these phase errors, but is limited 
by the timescale on which phase can be accurately measured on the target source itself. 

In the case of a Boo, the measured peak flux density in the normally calibrated image is about 
50 mJy. Given the system temperature (350 K), bandwidth (670 MHz, reduced from 800 MHz due 
to edge-channel flagging in the wideband average), and a correlator efHciency of 0.88 (2-bit) for 
these observations, the expected noise in a single-interferometer measurement is about 70 mJy in 
10 min. Using self-calibration with a 9-element array improves the SNR for antenna-based phases 
by about a factor of 3, due to baseline averaging. The expected r.m.s. noise in the antenna-based 
phases using self-calibration is, therefore, about 25 mJy in 10 min. Assuming the true flux of 
a Boo to be 80 mJy, we would expect an SNR of about 3 using 10-min averaging on the source. 
The expected phase error for an SNR of 3 is ~ 20° for each antenna (arcsin(l/3)). The decorrelation 
expected for an r.m.s. phase error of 20° is about 6% (reduced amplitude = exp(— r^/2), where 
r=the r.m.s. phase error in radians). It would, therefore, appear safe to use an interval of 10 min 
or longer to self-calibrate the 1-mm a Boo data in phase. Of course, any additional phase error 
(e.g., atmospheric) incurred during the self-calibration interval will increase the decorrelation in 
the final image. 

We can estimate the amount of atmospheric decorrelation expected. The average r.m.s. delay 
path in 10 min for a 100 m baseline on the day of the a Boo 1-mm observations (01 Jun 2000) 
was 210 /xm, or ~1 rad at 1360 /xm wavelength. In D-array, the maximum baseline is about 30 m. 
Assuming that the atmospheric r.m.s. phase scales as (baseline) the r.m.s. phase becomes 
0.37 rad on the longest baseline [(30 m/100 * 1 rad], and about 0.2 rad for the median baseline 

of 15 m. The expected atmospheric decorrelation in 10 min for the 30 m baseline (r=0.37 rad) is 
0.85, and for the 15 m baseline (r=0.2 rad) is 0.98. Thus, atmospheric decorrelation on 10-min 
timescales is not expected to be very important, and phase variations on the longer timescales will 
be removed by the self-calibration. 

We validated these ideas on the June 2000 observations of a Boo by experimenting with self- 
calibration timescales in the reduction of the lower side-band data. The total cleaned flux in the 
resulting images of this star showed a spread of only 3.5%, while the more relevant quantity, the 
stellar peak flux densities, were 82.1, 81.5, and 81.1 mJy for a 15-, 20-, and 30-min timescale. We 
settled on a 20-min self-calibration timescale to reduce the data for a Boo, and 30 min for a Tau. 

The final rows in Table 1 present the combined peak flux density information based on the 
totality of the available, good quality, stellar data for a Tau and a Boo. The corresponding 
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combined images from which these final flux densities were determined appear as Figure 1. 

4. Observations with LWS on ISO 

Observations of a Tau with full wavelength coverage were made in three revolutions: 818 
(TDT 81803302), 848 (TDT 84801701), and 861 (TDT 86101101). The observations of a Tau each 
consisted of 28 scans, 14 in each direction with a spectral oversampling of 4. Each observation 
had about 4630 s integration. LWS made two AOT observations of a Boo in revolutions 448 
(TDT 44800303) and 608 (TDT 60800904), each with about 2300 s integration. These observations 
consisted of 12 scans, 6 in each direction, taken with an oversampling of 4. 

Each LWS spectrum was processed with Offline Processing version 10 (OLPIO) and all the data 
were reduced using the Infrared Spectral Analysis Package (ISAP). The reduction steps consisted of 
removing outliers and averaging over individual detectors at maximum resolution, clipping at 2.5(7. 
No attempt was made to scale the detectors to make the sub-spectra overlap in the reductions. We 
accomplished this outside ISAP using the splicing techniques described by Cohen et al. (1992b) for 
constructing absolutely calibrated spectra. Care was taken to use the same number of scans in each 
direction, so the effects due to transients in the data are minimal, and no transient correction was 
applied. The fundamental Uranus model for calibrating the LWS is a synthesis of globally averaged 
Voyager InfraRed Interferometer Spectrometer (IRIS) data up to 50 /im and James Clerk Maxwell 
Telescope data covering the 0.35-2 mm wavelength range (Griffin & Orton 1993) with a radiative 
transfer model similar to the one described by Griffin & Orton (1993). 

A formal off-source, full-grating, LWS integration is available only for a Boo (TDT 60801005), 
and two fixed-grating (sampling only ten wavelengths with no detector scans) observations exist 
(TDTs 21701113 and 27202302). However, it is vital to recognize that corrections for off-source 
emission in LWS are inappropriate for faint objects like normal stars in low-cirrus regions. 

All observed LWS spectra have had the instrumental dark currents subtracted from them 
during the reductions. In our OLPIO reductions, we used the "fixed dark currents", as determined 
in dedicated calibration observations^ . These LWS dark current measurements were essentially dark 
background measurements. The Fabry-Perot was placed in the beam and offset. Therefore, there 
was flux on the detectors, radiated by the 3.5-K blank, but this was negligible. The signals measured 
remained constant throughout the mission, and were subtracted from the measured photocurrents 
on-source. When the dark signal, measured from the blank, was compared with the signal measured 
on the dark patch of sky, no difference was found. Therefore, for dark sky, no additional "sky" 
subtraction is necessary. For regions with high cirrus, there is a sky signal in the LWS beam and, if 
a point source is being observed, an "off' position is needed to determine the local sky background. 
In dark regions, the sky flux predicted by COBE can form a substantial fraction of the LWS dark 



^http://www.iso. vilspa.esa.es/manuals/HANDBOOK/lws_hb/node30. html 
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signal. However, this is deemed undetected, because no additional signal is detected over the signal 
from the blank. 

Consequently, we must assess the measured sky background around the two K giants at those 
times when the LWS spectra were taken. If those backgrounds were below the fixed dark values, 
then there is no need for further subtraction of off-source LWS spectra. Indeed, that subtraction 
would produce meaningless stellar spectra by removing the dark currents twice. Both LWS and 
ISOPHOT, off-source, measure the sum of zodiacal (hereafter "zodi"), cirrus, and extragalactic 
emission, as did DIRBE. "IRSKY"^ shows that extragalactic and cirrus confusion noises are far 
below zodi, so the dominant background for the two K giants is zodi. We have used two methods 
to evaluate the FIR background near a Tau and a Boo: actual measurements in off-source pixels 
when these stars were observed in the FIR arrays of ISOPHOT; and, as a coarse check, DIRBE 
predictions for the LWS observations based on their solar elongation angles. 

The computation of ISOPHOT background sky measurements near a Tau and a Boo is sum- 
marized in Tables 2 and 3. These present the relevant ISOPHOT observation by filter, nominal 
wavelength, TDT and measurement number, and the background brightness and its 1-cr uncer- 
tainty. Our conversion from point source flux to surface brightness is different from the one given 
in the ISOPHOT handbook (Laureijs et al. 2003). The supporting analysis and justification are 
given by Schulz^. When multiple data exist for a filter, inverse- variance weighted combinations of 
those data arc also given. Conversion from MJy sr^"*^ to Fa (in W cm~^ ^m~'^) was accomplished 
using the published values for LWS detector beam solid angles, effective apertures, and bandwidths. 
Table 4 presents these backgrounds, and compares them with the fixed dark currents and DIRBE 
predictions (interpolated between 140 and 240 fim to assess the 160- fim values), matching each 
ISOPHOT band to the nearest LWS detector. The DIRBE predictions are generally in accord 
with ISOPHOT measurements. Clearly, the measured sky backgrounds around the two stars are 
all below the fixed dark currents, confirming that only low level cirrus appears in their vicinities 
(with the sole exception of a Boo in LW2, where ISOPHOT suggests 50% more sky brightness 
than the LW2 dark current). The DIRBE predictions are generally in accord with the ISOPHOT 
measurements. Subtraction of off-source LWS spectra would be inappropriate for these K giants. 

Figures 2 and 3 illustrate our reduced, spliced LWS spectrum for each star, presenting the 
mean spectrum over all TDTs, with the ibScr bounds for a Tau (which has more LWS data, longer 
integrations, and higher SNR) and ±2(7 bounds for a Boo (with its noisier data sets). The sizeable 

triangular excursions in both spectra near 53 and 106 /xm represent the difficulty, for NIR-bright 
normal stars, of mitigating the substantial leaks of 1.6-/im radiation to which LWS is prone. We 
have cleaned the observed spectra by 30-point boxcar smoothing of the data (equivalent to about a 
2-|Ltm interval), a Tau's combined spectrum contained 1181 wavelength points; a Boo's 706 points. 



^http://www.ipac.caltech.edu/ipac/services/irsky/irsky.html 
^http://www.iso. vilspa.esa.es/users/expLlib/PHT_list.html/fpsf_report.ps.gz 
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Each figure incorporates the absolutely calibrated model continuum spectra that we have computed 
(sec below) as the almost flat long-dashed lines that cross the figures, and the meanzblcr bounds 
on our older calibrated energy distributions (dash-triple-dotted lines) delivered in January 1996 
for the calibration of ISOPHOT. To convert the model surface fluxes to predicted fluxes at the 
Earth, we have adopted the angular diameter, 20.88=b0.10 mas, of Ridgway et al. (1982) (using 
lunar occultation) for a Tau and the angular diameter, 21.0±0.2 mas, of Quirrenbach et al. (1996) 
for a Boo (measured by optical intensity interferometry). 

Our new, more detailed, calculations of stellar continuum spectra fall well within the uncer- 
tainties associated with the calibrated spectra originally delivered to ISO. We wish to quantify the 
degree of accord between measured and modeled spectra in the FIR for each star. The grating 
simultaneously illuminated the ten LWS detectors, and each detector was separately wavelength 
and flux calibrated. Therefore, LWS consists of ten spectrally independent detectors. The degree 
of oversampling of these stellar spectra was 4; consequently, adjacent points in Figures 2 and 3 
are not independent. Therefore, to assess the consistency of observations and predictions, and to 
handle the dependence of adjacent LWS points in the fully-plotted spectra correctly, we include the 
data binned by detector as a series of bold crosses in Figures 2 and 3. Filled circles represent the 
inverse-variance weighted average of all the meaningful data in each LWS detector (e.g. after the 
rejection of data contaminated by if-band leaks). Error bars show the associated \-a uncertainties. 
Asymmetric horizontal error bars represent the actual wavelength ranges usable for each detector 
and star. Filled circles are plotted at the average wavelength of each detector's set of usable data 
for that star. 

The binned detector data lie within \a of the newly modeled energy distribution for a Tau, 
and within \o for a Boo. There is good agreement between the underpinning Uranus calibration 
of LWS, the totally independently calibrated stellar spectra made from the empirical spectra by 
our older model extrapolations based on the available observations between 8 and 23 /im, and the 
newly calculated, absolutely calibrated, emergent stellar energy distributions. Averaging the results 
for the independent detectors for each star with inverse- variance weighting, we determine that, for 
a Tau, this mean scale factor is l.OOOibO.Oll (over nine detectors) while, for a Boo, this factor is 
1.012ib0.025 (over seven detectors). Therefore, empirical stellar and planetary calibrations used by 
ISO, and stellar modeling, are absolutely reconciled to 3.3% (30") for the better measured a. Tau, 
and 7.5% (3a") for the more poorly measured a Boo. 



5. The Radiative Equilibrium Model Atmospheres 

5.1. Atmospheric parameters and computation of the emergent spectrum 

The spectrum calculations were carried out with the SOURCE model atmosphere program in 
essentially the same fashion as described in Carbon et al. (1982). Speciflcs pertinent to the current 
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calculations are described in this section. 

For a Tau, we adopted an effective temperature of 3920 K (Blackwell et al. 1991) and a 
log(g) of 1.5 (Smith k. Lambert 1990). For the isotopic abundances of C, N and O, we adopted 
log e(C) = 8.40, log e(N) = 8.20, log e(0) = 8.78, i^c/i^c = 10, ^^O/^^O = 562, and ^^O/^^O = 475 
(Harris & Lambert 1984; Smith & Lambert 1990). The remaining metal abundances were chosen 
to be solar, following the lead of Smith k. Lambert (1985). The preceding abundances were de- 
rived using the MOOG spectrum synthesis code (Sneden 1973) and model atmospheres (Johnson 
et al. 1980) constructed with the ATLAS (e.g., Kurucz (1993)) model atmosphere code and opacity 
sampling. Our choices of Tg//, log(g), and abundances are entirely consistent with current deter- 
minations of the parameters of a Tau. In particular, the values are well within the error bars of 
the values recently derived directly by Decin and her collaborators (Decin et al. 2000, 2003) from 
fits to stellar spectral energy distributions generated by using the MARCS (Gustafsson et al. 1975; 
Plez et al. 1992) model atmosphere code. Bertrand Plez (Plez 1999) kindly computed an a Tau 
model atmosphere for us using the aforementioned parameters and the programs described in Plez 
et al. (1992). Thus, this model was computed with essentially the same model atmosphere program 
and opacities as employed by Decin. 

For a Boo, we adopted a slightly revised version (Peterson (1995): hereafter PDK) of the 
Peterson et al. (1993) model atmosphere which has more depths than the original published model. 
This atmosphere has an effective temperature of 4300 K, log(g) of 1.5, [Fe/H] = -0.5, and other 
elemental abundances as determined by Peterson et al. (1993). The PDK model parameters derived 
by Peterson et al. (1993) were derived entirely from the observed stellar spectrum. As was the case 
for the Decin et al. investigation of a Tau, no assumptions were made regarding stellar temperature 
or gravity. We adopted all of the Peterson et al. (1993) values for our study. The parameters 
Tg//, gravity, microtiirbulent velocity, and [Fe/H] derived by Peterson et al. (1993) for a Boo are 
essentially identical to those derived later and independently by Decin et al. (2003), while the C, 
N, and O abundances agree within the error bounds. 

The alert reader will observe that we have adopted model atmospheres for a Tau and a Boo that 
were derived using entirely different model atmosphere codes (ATLAS and MARCS). We shortly 
will describe results obtained by using these atmospheres as input to a third model atmosphere 
code, our SOURCE. We do not feel that significant improvement would result from insisting that 
only ATLAS or MARCS models be used in our study. All three codes have been vetted against 
the solar case and should give very comparable results in the G-K star domain. We cite as support 
for this stance the striking agreement between the results of Peterson et al. (1993) and Decin et al. 
(2003) for a Boo. This occurs despite Decin et al. (2003) having used a completely different model 
atmosphere program and blanketing representation (MARCS with opacity sampling) to compute 
their a Boo model than used by Peterson et al. (1993) (ATLAS with opacity distribution functions). 
The agreement is even more impressive when one considers that these two investigations derived 
their results from two entirely different spectral regions, 0.5-0.9 ^um versus 2.38-12.0 /im. Wc regard 
this as evidence that different well-vetted model atmosphere codes and opacity representations 
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now produce comparable model atmospheres and visible-infrared spectra for K-giants within the 
limited approximations of static, homogeneous atmospheres, LTE, and mixing-length-theory 
radiativc/convective equilibrium. That docs not mean that the resultant models are accurate, or 
even adequate, representations of a real star. Differences between observed and computed spectra 
are now much more likely to arise from the presence of physics {e.g., NLTE, non-radiative heat 
deposition, and inhomogeneous dynamical structure) not treated in the simple models we use 
here. We will discuss this point in more detail as we proceed. 

The column mass (g cm"~^) vs. temperature relations for both the Plez a Tau and PDK a Boo 
models were linearly extrapolated to shallower depths as required to obtain optical transparency at 
the longest wavelengths. For the Plez a Tau model, the extrapolated layers are shallower than 0.020 
g cm~^; for the PDK a Boo model, shallower than 0.010 g cm~^. The resultant model atmospheres 

are shown in Figure 4. 

It should be stressed that both of the above adopted model atmospheres have purely radiative 
equilibrium surface temperature structures. Specifically, they do not have chromospheric temper- 
ature rises, a point whose importance will become evident later. We have included in Figure 4, for 
comparison, the Kelch et al. (1978) (KEL) a Tau chromospheric model B and the Ayres & Linsky 
(1975) (AL) a Boo chromospheric model. These chromospheric models are semiempirical, derived 
from non-LTE analyses of Ca II and Mg II line fluxes. While these chromospheric models are now 
over 20 years old, they are very similar from the upper photosphere, through the temperature min- 
imum region, and into the lower chromosphere to models of a Tau and a Boo used by Wiedemann 
et al. (1994), and the a Tau model derived by McMurry (1999). 

Using the Plez a Tau atmosphere and parameters detailed above, we have computed the 
synthetic spectrum for a Tau over the entire LWS range and beyond, 43 to 665 //m. The spectrum 
was computed using a full set of molecular line opacities which includes all the principal isotopomers 
of CO, SiO, OH, and H2O and using a variable wavelength mesh which locally provides f«2.5 points 
per Doppler width. To allow comparison with the LWS observations, we convolved the high- 
resolution result with a Gaussian instrumental profile function having a FWHM equal to 0.29 /im 
for A < 90.5 jim, and 0.60 /im, A > 90.5 /xm (e.g., Gry et al. (2003)). The resulting spectrum, 
normalized to the continuum and on a greatly magnified ordinate scale, is shown in Figure 5. A 
detailed description of the displayed spectrum may be found in Carbon et al. (2005). We do note 
here that almost all the absorption beyond 50 fim is due to the isotopomers of SiO. At shorter 
wavelengths the principal contributor is OH with some H2O and CO as well. 

It is clear that, at the ISO LWS and the BIMA spectral resolutions, all differences between the 
pure continuum and the synthetic spectrum computed with the full set of molecular line opacities 
are less than 2%; beyond 130 nm, they are less than 1%. This is much less than the noise level 
in the observed a Tau spectra in this region. As a consequence, there is no significant advantage 
gained by using the full synthetic spectrum. For clarity in our presentation, we will use only the 
continuous spectrum in our comparisons with the a Tau observations. Since a Boo is both hotter 
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and has lower metallicity, we shall use the continuum alone for comparisons there as well. 

5.2. Depth of formation and brightness temperatures 

For both the Plez a Tau and the PDK a Boo models, the dominant continuous opacity in the 
continuum- forming layers for wavelengths redward of the H~ opacity minimum at 1.64 ^m arises 
from H~ free-free. This opacity increases in strength as A^. As a result, the emergent continuum 
flux is formed at progressively shallower layers in the stellar atmosphere as one proceeds to longer 
wavelengths. Figure 4 illustrates this point for the Plez a Tau and the PDK a Boo models. In this 
figure we have marked, for the radiative equilibrium models, the atmospheric depths for which the 
continuum optical depth, r|, reaches unity at the indicated wavelengths. That is, we have marked 
the column masses at which = 1.0; roughly speaking, in the FIR, depths near this value and 
shallower make the greatest contribution to the emergent continuum flux at wavelength A. Note, for 
example, that the continuum fluxes between 1 mm and 3 mm are sampling very different parts of the 
atmospheric temperature structure in both radiative equilibrium models than, say, the continuum 
fluxes at 2.3 /xm. Indeed, the longest wavelengths are sampling, in the continuum, atmospheric 
levels reached only at the cores of lines in bluer portions of the spectrum. Only when one proceeds 
to the vacuum ultraviolet and its strong photoionization continua of metals is it possible to reach 
in the continuum the shallow atmospheric levels probed by the continuum at long wavelengths (e.g. 
Figure 1 of Vernazza et al. (1976)). Thus, the FIR continuum is an effective probe of the shallowest 
stellar layers which can supplement information derived from line cores and UV continua. 

As noted above, there is very little difference between the full synthetic spectrum and the 
spectrum of the continuum. Figure 4 provides an important clue as to why spectral lines generally 
become less prominent at longer wavelengths. It is evident that, as one goes to shallower depths in 
the models, the slope of the temperature-column mass relation becomes smaller. The significance of 
this can be appreciated by considering an absorption line in the a Tau model with a characteristic, 
depth-independent, line core/continuum opacity ratio of 1000. If the absorption line occurs at a 
wavelength of 2.3 //m, the continuum will be formed at ~5200 K and the line core will be formed at 
?a3200 K. The same line at a wavelength of 200 /xm will have the continuum formed at Ri3100 K and 
the line core at !=a2600 K, a much smaller temperature difference between line core and continuum 
than at 2.3 //m. Couple this with the decreasing temperature sensitivity of the Planck function 
as one goes to longer wavelengths and one has a recipe for progressively less prominent spectral 
features. 

Figure 6 provides another window into the wavelength dependence of continuum formation. In 
the upper panel of Figure 6 we show the column masses for which = 1.0 over the spectral range 
0.5 nm to 200 //m. As noted above, there is a steady march to progressively shallower depths of 
continuum formation as one goes to longer wavelengths. Since the dominant H" bound-free opacity 
peaks at roughly 0.825 /xm, continuum fluxes at wavelengths between 0.825 /xm and 1.64 /xm probe 
the same atmospheric layers as continuum fluxes between 1.64 /xm and 6.5 /xm. Longward of 
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6.5 fim, the continuum flux probes atmospheric layers accessible only in line cores and in the metal 
bound-free continua of the vacuum ultraviolet. 

In the lower panel of Figure 6 we show the continuum "brightness temperature" over the 
same spectral range for both models. The brightness temperature is defined as the temperature 
of a blackbody which gives the same flux as the model atmosphere at the indicated wavelength 
(a precise definition is given in Section 6.3). Three points are important. First, neither a Tau 
nor a Boo could have its continuum fluxes over any large spectral range matched by a blackbody 
of a single temperature. Second, the blackbody temperatures that do match the model fluxes at 
speciflc wavelengths are generally quite different than the effective temperatures of the models (or 
the stars) themselves. Third, comparing with Figure 4, one sees that the brightness temperature is 
systematically smaller than T(r^ = 1.0), the temperature of the atmospheric layer at which optical 
depth unity is reached in the continuum at wavelength A. (For an atmosphere with a linear source 
function at wavelength A, the brightness temperature at A and T(r^ = 2/3) would be identical. 
Unfortunately, in the particular cases of the PDK and Plez models, the FIR source functions 
are not linear in but, rather, quadratic.) The differences between T(r^ = 1.0) and brightness 
temperature vary systematically with wavelength. The differences are greatest at 1.64 /im and 
decrease to the red. Beyond 40 /um, the differences are always less than 100 K. (Had we chosen 
for illustration the more correct, but still approximate, T(t^ = 2/3), the temperature differences 
would be smaller yet.) Thus, the brightness temperature is a good flrst guess at T(r^ k. 1.0) at the 
longer wavelengths. Later, when we discuss differences between observed and predicted brightness 
temperatures, the reader may use this approximation and Figure 4 to roughly estimate the implied 
errors in the radiative equilibrium model temperature structures. 



6. Discussion 

6.1. Evidence for significant fiux excesses at 1.4 mm and 2.8 mm 

Table 5 shows the comparison between the observed fluxes and the predictions of the radiative 
equilibrium models. Only for the case of a Tau at 1.4 mm does there appear to be approximate 
agreement between the observed flux and the flux predicted from the radiative equilibrium model. 
The remainder of the observed fluxes are all significantly larger than the fluxes predicted by the 
models. The discrepancies are particularly striking in the case of a Boo. 

We do not believe that the fiux computation for the radiative equilibrium models could have 
errors sufficiently large to produce the marked differences between computed and observed mil- 
limeter fiuxes. Our model atmosphere code has been carefully vetted and the dominant continuous 
opacities are themselves quite well determined (e.g., van der Bliek, Gustafsson, &; Eriksson (1996)). 
(As an example of evidence for this, we refer the reader forward in the text to Figure 11. There we 
show the continuum for the Plez a Tau model computed with an entirely independent model code.) 
On the observational side, our derived errors for the observed 1.4-mm and 2.8-mm fluxes preclude 
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dismissing the differences between observed and computed fluxes as simply due to measurement 
error. We believe other explanations must be considered. 

Since the adopted angular diameters for these two stars appear to be extremely well determined, 
with standard deviations of ~1% or less, they would not seem to be a significant source of error. 
Moreover, a major change in the angular diameters would compromise flux agreement for the A < 
200 iim region. Nevertheless, it should be noted that the Quirrenbach et al. (1996) diameter is 
an arithmetic mean of values measured from 0.45 /xm to 2.2 /zm and the Ridgway et al. (1982) 
diameter is a error- weighted mean biased toward values obtained at 1.6 /xm and 3.8 fxm. At 
any particular monochromatic wavelength, an observed angular diameter represents the apparent 
diameter of the ~ 1.0 surface of the star. Since, as shown in Figures 4 and 6, there is a progressive 
march in depth of continuum formation to ever shallower atmospheric layers as one moves redward 
of 1.64 iim, the angular diameters themselves should increase with wavelength. If the angular 
diameters at 1.4 mm and 2.8 mm are substantially larger than those at the wavelengths of the 
Quirrenbach et al. (1996) and Ridgway et al. (1982) observations, our predicted fluxes in Table 5 
could be seriously underestimated. We can roughly estimate the magnitude of this source of error, 
at least in the case of the radiative equilibrium models. In Tabic 6 wc show the physical radii of 
a Tau and a Boo based on the indicated angular diameters and Hipparcos parallaxes (ESA 1997). 
We also show the physical distances in the radiative equilibrium models between rf 64^^ ~ ^-^ 
t"^ —in 

'2.8mm ~ 

Given the tabulated numbers, it is clear that, for the PDK and Plez models, the atmospheric 
extent in the continuum at 2.8 mm is not sufficient to seriously affect our predicted fluxes. We 
find a similar result in the case of the only model chromosphere with a published height scale, 
McMurry's chromospheric model for a Tau (McMurry 1999). The McMurry atmosphere is plane- 
parallel and hydrostatic, constructed with log g = 1.25. Note that this gravity leads to a more 
distended atmosphere than docs the higher gravity we adopted. The McMurry atmosphere has 
an extent of ~2.6 x 10^^ cm between the layers where the 1.64 jim continuum is formed and the 
layer with T = 7120 K in the chromosphere. This extent is only 8% of the stellar radius and again 
seems too small to explain the discrepancy between computed and observed fluxes. Finally, as 
we will describe in Section 6.3, Drake Sz Linsky (1986) determined that a Boo has essentially the 
photosphcric radius at 2 cm. Given the expected increase in angular diameter with wavelength, 
this would imply that a Boo would have the photospheric radius at < 2 cm as well. 

Based on the above arguments, we believe that the flux excesses observed at 1.4 mm and 
2.8 mm, relative to the predictions of the purely radiative models, are real and significant. Fur- 
thermore, as we shall show, there is good evidence that these excesses are signatures of the chro- 
mospheric temperature rise in the two red giants. 
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6.2. Other observations of the sub-mm through cm region 

Figures 7 and 8 present our millimeter flux densities (filled circles) and the color-corrected 
IRAS measurements of each star at 12, 25, 60, and 100 //m (crosses). Each figure includes the Plez 
a Tau or PDK a Boo spectrum, as appropriate, from 10 /xm to 3 mm. 

There are additional high-frequency radio continuum observations of these two K giants in the 
literature. Altenhoff et al. (1986) made single-dish 86-GHz measurements of a Boo (21.4±7.5 mJy) 
with the IRAM 30-m telescope, while Altenhoff et al. (1994) likewise detected both a Tau and 
a Boo at 250 GHz, obtaining 51ib6 and 78ib8 mJy, respectively. There are also VLA centimeter 
measurements on both stars at 4.9, 8.4, and 14.9 GHz by Drake & Linsky (1986), and unpublished 
data by Drake, Linsky, & Judge from the compilation of radio data by Wendker (1995). These 
long wavelength measurements all have substantial uncertainties, in part due to the faintness of the 
stars, in part to the potential for confusion with extragalactic background objects. Nonetheless, in 
combination with our own radio data (Figures 7, 8) one sees that these high-frequency data (open 
diamonds) are in good agreement with our aperture synthesis measurements (filled circles). 



6.3. Brightness temperatures for all the observed fluxes 

It is difficult to fully appreciate the significance of the excess fluxes from the logarithmic 
scales of Figures 7 and 8. An insightful alternative view may be obtained by combining the mm 
and cm fluxes from other observers with our data and looking at the whole dataset as brightness 
temperatures over a broad spectral range. In Figures 9 and 10 we show the brightness temperatures 
corresponding to all the observed flux data between 43 ^um and 6 cm used in our previous figures. 
For comparison we have also plotted the predictions of the radiative cqTiilibrium models. In deriving 
brightness temperatures from the observed fluxes, we used the angular diameters of Table 6 in the 
following relation: 

^ 14387.75 /A 

where Sjy{\) is the observed flux in Jy, a is the angular diameter in mas, A is the wavelength 
in microns, and Tb(A) is in Kelvins. The error bars in Figures 9 and 10 correspond only to the 
uncertainties in the observed fluxes. In the case of the a Tau data at 2 cm and 3.6 cm, no error 
values were reported by Wendker (1995) and we used our own rough estimates of the plausible 
errors. Before proceeding, a caveat needs to be noted regarding the longer wavelength points 
especially. The observed fluxes presumably come from progressively shallower stellar layers as A 
increases. All of the brightness temperatures in Figures 9 and 10 have been deduced assuming 
the photospheric angular diameters of Table 6. Atmospheric extension may play a role in the 
chromospheric layers of a Tau and a Boo. Carpenter et al. (1985) estimated the extension of the 
chromospheres of these stars. They argued that the C II UV 0.01 lines which they studied were 
formed in chromospheric layers with temperatures of 10000 K in a Tau and «8000 K in a Boo. 
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They deduced that these chromospheric layers had radii of f«2.0 R* and Ril.4 R*, respectively, 
where R* is the stellar photospheric radius. Although these estimates are, in the authors' words, 
"relatively crude", they do suggest that atmospheric extension may play a role at least in the 
hotter layers of the chromosphere. Different evidence of extension was found by Drake &: Linsky 
(1986) in their study of 2 cm and 6 cm fluxes. They treated all the radio emission as coming from 
an optically thick ionized wind for which a "half-power radius" could be deduced. In the case of 
a Boo they found the 2 cm half-power radius was paR^; at 6 cm, the half-power radius had grown 
to ~1.7 R*. Taking the long wavelength approximation to the brightness temperature expression 
above, we see that, for a given observed flux, Tb(A) oc where a is the angular diameter. Since 

the shallowest layers may have angular diameters significantly larger than those of the photosphere, 
the corresponding brightness temperatures in Figures 9 and 10 should be regarded as upper limits. 

6.4. Evidence for contributions from a stellar chromosphere 

We now describe why we believe that the flux excesses at the longer wavelengths arise from 
the presence of stellar chromospheres in a Boo and a Tau. 

6.Jt..l. Circumstantial Evidence 

Circumstantial evidence suggests the possibility that stellar chromospheres could explain the 
flux excesses. This evidence comes by examining the depths at which = 1.0. In the case of a Boo, 
as shown in Figure 4, the AL chromospheric model departs from the purely radiative equilibrium 
model at a column mass of ~ 0.75 g cm^^ (Ayrcs &: Linsky 1975). This is just in the column mass 
range of the PDK radiatve equilibrium model where optical depth unity is reached between 200 /um 
and 1 mm. Given this, it is not surprising that the PDK model might significantly underpredict 
the flux at both 1.4 mm and 2.8 mm in a Boo. We believe that the progressively larger differences 
between the predictions of the PDK model and the observations beyond 120 fj,m are due to the 
contributions of a Boo's chromosphere. 

The circumstantial evidence in the case of a Tau is not so clear. As shown in Figure 4, the KEL 
a Tau chromospheric model begins its temperature rise at a column mass of ~ 0.3 g cm^^ (Kelch 
et al. 1978). This is slightly shallower in the Plez a Tau model than where mm ~ ^-^ occurs. 
Nevertheless, the 2.8 mm observed brightness temperature is distinctly higher than that predicted 
by the Plez model in Figure 9. In fact, with the exception of our 1.4 mm point, the general run 
of observed points beyond 200 /um suggests that some enhanced (chromospheric?) emission may 
be present at 1 mm or even shorter wavelengths. Whether this is true or not will take additional 
observations in this spectral region with better SNR. 
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6.4-2. Evidence from NLTE modeling 

While NLTE chromospheric modeling is well beyond the intent of this investigation, one of us 
(ADM) has computed the infrared fluxes predicted by his chromospheric model of a Tau . Details of 
the model and the NLTE computation may be found in McMurry (1999). This model is essentially 

identical to the KEL model from the photosphere through the lower and middle chromosphere, 
differing from the KEL model only at column masses shallower than 10~'^ g cm^^ (sec Fig. 5, 
McMurry (1999)). The NLTE radiative transfer computation was carried out with the independent 
code MULTI. We show in Figure 11 the results of this NLTE computation with a chromospheric 
model compared with the observations, and with the predictions of the Plez model. 

It is immediately apparent that, again with the exception of our 1.4 mm point, the observed 
data beyond 200 /xm is better fit by the NLTE chromosphere model than the LTE Plez radiative 
equilibrium model. We believe this is good evidence that the stellar chromosphere of a Tau has a 
significant impact on the observed fluxes beyond 200 /xm. Indeed, it should be noted that the fit is 
even improved between 100 /xm and 200 /xm. 

6.5. The transition from photosphere to chromosphere 

For a Boo (Figure 10) it appears that the transition from purely radiative occurs beyond 
120 ^m, and certainly well before 1.4 mm. In our opinion this makes any purely radiative model 
for a Boo a poor choice for calibration in the sub-mm spectral region. This appears to be the 
case for a Tau as well, especially considering the predictions of the McMurry chromospheric model 
(Figure 11). Given the evidence described above, we do not believe that it is prudent to use cool star 
radiative models for calibration purposes at the longer wavelengths unless alternatively calibrated 
observations already confirm the validity of the purely radiative models. One must be certain 
that the flux predicted by the purely radiative model at the long wavelength end of the region 
being calibrated by the theoretical stellar spectrum does indeed correspond to the actual star. In 
other words, we believe that when one is approaching the spectral regions where a chromospheric 
contribution can be significant, one can safely use purely radiative cool star models only as an 
interpolation tool for calibrations, not as an extrapolation tool. 

We would like to emphasize the importance of the spectral region from psSO /xm to several 
millimeters to furthering an understanding of cool stellar atmospheres. Fluxes from this spectral 
region reflect the mean time- and spatial-averaged emission from the upper photosphere and chro- 
mosphere for stars like those we consider here (see Figure 4). Accurate fluxes measured through this 
spectral region should allow one to test directly model atmosphere structures, even complex ones 
that include NLTE and dynamics, over the very interesting depth range in which non-radiative 
processes like shocks begin to alter the structure. By selectively sampling progressively longer 
wavelengths, one can follow the important transition from photosphere to chromosphere, locating 
the depth at which purely radiative models like the PDK and Plez models begin to fail. 



-18- 



6.6. Stellar absolute calibrators for ISOPHOT 

Our original, theoretically extrapolated, continuum spectra for the two K giants, supplied to 
ISOPHOT as absolutely calibrated energy distributions, are in extremely good agreement with our 
new calibrated synthetic spectra (Figures 2, 3), well within the ±6% uncertainties assigned to the 
original calibrators (Cohen et al. 1996b). The LWS spectra, reduced with respect to Uranus as 
calibrator, are consistent with the new (and old) synthetic spectra, thereby uniting FIR planetary 
and stellar calibrations at the few percent level (Figures 2, 3). 

There are independent data on a Boo, derived from the validation of ISOPHOT's final pipeline 
(OLPIO) (Richards 2003), that show it conforms to our 1996 model spectrum. In addition, Schulz 
(2003) has described the performance through the ISO mission of a significant number of ISOPHOT 
calibrators. By removing the contributions made by any particular source to the calibration process, 
and calibrating that source with respect to all other calibrators, one can assess the ratio of observed 
to modeled irradiance. Using the database described by Schulz (2003), we have been able to extend 
the range over which our new model of a Boo can be tested, and to validate q Tan independently 
of our LWS data. For a Boo, we have determined the observed to modeled irradiance ratios for 
the ISOPHOT ClOO, C120, C135, C160, C180, and C200 filters ("C" denotes the FIR camera 
within the ISOPHOT instrument, and each filter is designated by its reference wavelength in \xm). 
Expressing these observed irradiance ratios with respect to our new modeled spectrum for a Boo 
we find: ClOO, 1.02±0.02; C120, 1.14±0.08; C135, 1.01±0.02; C160, 1.06±0.02; C180, 1.00±0.03; 
C200, 1.00=b0.26. For a Tau, measurements are available solely in the C160 filter and yield a ratio 
of 0.96ib0.05, with respect to our new model. We conclude that the separate absolute calibrations of 
LWS and ISOPHOT are in accord within the uncertainty of the observations, despite their different 
dependences on planets, asteroids, and stars. This should not be interpreted to mean that more 
accurate observations might not reveal discrepancies between FIR observations and the predictions 
of static, LTE model atmospheres computed in simple radiative/convective equilibrium. 



7. Conclusions 

Our primary objective was to investigate whether purely radiative models for two normal K 
giants could reproduce the observed FIR and microwave fluxes. We have found that for both a Tau 
and a Boo the radiative equilibrium temperature structures of the upper photospheres produce 
fluxes in agreement with the observations, but only up to a certain point. For a Boo and a Tau the 
agreement extends to at least 125 /xm, corresponding to a column mass of Ri5 gm cm~^ in the PDK 
atmosphere and ?»10 gm cm~^ in the Plez atmosphere. Both K giants clearly show definite chro- 
mospheric emission at 2.8 mm. In a Boo this is evident by 1.4 mm. Only additional observations 
can pinpoint precisely where between 125 /im and 2.8 mm the chromospheric contribution becomes 
significant. For this reason, we cannot rely on purely radiative model atmospheres to form a basis 
for stellar calibration in the sub-mm region. Purely radiative models of the continuum might be 
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used to provide calibrations between those shorter FIR wavelengths where the validity of the model 
fluxes already has been independently verified by observations. 

Finally, our results indicate that accurate FIR, sub-mm, and mm observations can provide 
valuable new information for researchers modeling stellar chromospheres of cool stars. 
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Fig. 1. — Images of both stars at 3 mm and 1 mm constructed from the combination of all available 
data. In each plot, contours are given at -20, -10, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100% of the 
peak flux density, a) a Tau at 3 mm with a 13.97 mJy peak, b) a Tau at 1 mm with a 25.78 mJy 
peak, c) a Boo at 3 mm with a 20.09 mJy peak. The offset of the peak from the stellar position 
after allowance for proper motion is caused by the process of self-calibration, d) a Boo at 1 mm 
with a 83.5 mJy peak. 
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Fig. 2. — LWS spectrum of a Tau showing the mean observed spectrum (sohd line) after boxcar 
smoothing, and the meanztSa bounds (short-dashed hues), compared with the Plez model con- 
tinuum spectrum (long-dashed lines) and with the meaniblcr bounds on the old SED provided to 
ISOPHOT in 1996 (dash-triple-dotted lines). Each detector result is shown as a heavy cross, where 
horizontal bars represent the actual wavelength extent of the useful data, and vertical bars show 
the l-o" errors on the plotted weighted-mean points. 
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Fig. 3. — LWS spectrum of a Boo showing the mean observed spectrum (soHd Hne) after boxcar 
smoothing, and the nicanib2(T bounds (short-dashed hues), compared with the PDK model con- 
tinuum spectrum (long-dashed lines) and with the meanilcr bounds on the old SED provided to 
ISOPHOT in 1996 (dash-triple-dotted lines). Heavy crosses as for Fig. 2. 
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Fig. 4. — The adopted model atmosphere temperature structures for a Boo and a Tau. Arrows 
mark the depths at which optical depth unity is reached in the continuum at the indicated wave- 
lengths. The upper set of arrows refers to the PDK model; the lower set, to the Plez model. For 
comparison, the AL and KEL chromosphere models are shown as dashed lines. 
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Fig. 5. — Continuum normalized flux predicted by the Plez a Tau model for the spectral range 
43 //m to 665 jim.. Note that the vertical scale is greatly magnified. 
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Fig. 6. — Upper panel: The column masses at which continuum optical depth unity is reached 
in the PDK a Boo and Plez a Tau radiative equilibrium models. Lower panel: The brightness 
temperatures predicted by the PDK a Boo and Plez a Tau radiative equilibrium models. 
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Fig. 7. — FIR to cm spectrum of a Tau. Color-corrected IRAS points (crosses) are plotted together 
with our mm observations (filled circles), the LWS spectrum (solid), and our computed continuum 
spectrum (dashed line). Also shown are mm and cm data from the literature (open diamonds). 
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Fig. 8.— 



FIR to mm spectrum of a Boo as for Figure 7. 
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Fig. 9. — Brightness temperatures corresponding to the observed and computed fluxes for a Tau 
which have been discussed in this paper. LWS observations are represented by open squares, our 
mm measurements by asterisks, all other observations by open circles. Error bars represent l-cr 
deviations in the fluxes only; errors for 2-cm and 3.6-cm points are our own estimates (see text). 
Brightness temperatures predicted by the Plez model are shown by a solid line. Note that the 
brightness temperatTires at cm wavelengths may be only upper limits. See the text for a discussion 
of the impact of wavelength dependent angular diameters. 
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Fig. 10. — Brightness temperatures corresponding to the observed and computed fluxes for a Boo 
which have been discussed in this paper. LWS observations are represented by open squares, our 
mm measurements by asterisks, all other observations by open circles. Error bars represent \-g 
deviations in the fluxes only. Brightness temperatures predicted by the PDK model arc shown by 
a solid line. Note that the brightness temperatures at cm wavelengths may be only upper limits. 
See the text for a discussion of the impact of wavelength dependent angular diameters. 
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Fig. 11. — Brightness temperatures for a Tau including the predictions from the McMurry NLTE 
chromospheric model and from an independent computation using the LTE Plez model. Observa- 
tions are marked exactly as in Figure 9 with the same error bars. We have magnified the vertical 
scale slightly to reduce the clutter of points. The brightness temperatures predicted by the Mc- 
Murry NLTE chromospheric model are indicated as a dashed line. Brightness temperatures for the 
Plez model in LTE computed using MULTI are shown as open circles without error bars. Brightness 
temperatures predicted by SOURCE from the Plez model are shown by a solid line as before. Note 
that the brightness temperatures at cm wavelengths may be only upper limits. See the text for a 
discussion of the impact of wavelength dependent angular diameters. 
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Table 1. Journal of observations 



Star 


Date 


Array 


Frequency 


Time 


Phase 


Secondary 


Primary 


Traceable 


Peak F^ 








GHz 


on star 


Ref. 


Flux Ref. 


Flux Ref. 


to Primary 


mjy 


a Tau 


14JNovy7 




1 (JO. 0545 


4.32 hr 


0449+113 


3C111 


Mars, 
Uranus, Neptune 




17.14ibz.iy 


a Boo 


1 CAT 

loWov97 




1UD.D845 


4.22 hr 


1415+100 

1 Q ^7 1 1 CV) 

looY+iyo 


3C111 


Mars, 
Uranus, Neptune 




20.27+2.95 


a Boo 


15Nov97 


C 


110.1282 


4.22hr 


1415+133 
1357+193 


3C111 


Mars, 
Uranus, Neptune 




18.47+2.73 


a Boo 


15Nov97 


C 


Both 












19.73+2.16 


a Tau 


06Jun98 


c 


1UO.D845 


2.99hr 


0449+113 


3C111 


Maxs 




13.10+3.40 


a Boo 


llMar99 


B 


106.6845 




1357+193 


1415+133 




Mars 


24.57+1.20 


a Boo 


llMar99 


B 


110.1282 


4.37 hr 


1357+193 


1415+133 




Mars 


19.68+1.23 


a Boo 


llMar99 


B 


Both 












22.03+0.97 


a Boo 


12Mar99 


B 


10d.o845 


3.13 hr 


1357+193 


1415+133 


MWC 349, Mars 


Mars 


21.02+1.74 


cx Boo 


iziviaryy 


rS 


1 1 n 1 OQO 
liU. 1202 




loo ( +iyo 


1410+ioo 


iViWL^ o4y, Mars 


Mars 


lo.yUiti.zU 


a Boo 


izMaryy 


rJ 


riotn 












1 n no_l_i no 
iy.y2ihi.02 


CK Tau 


uooepyy 


■p» 
JJ 


2io.Uyc52 


i.ou nr 


UOoU+ioO 


UOoU+ioO 




iviwo o4y/iViars 


OA OQ-Lfi Q'? 

z4.zoitD.cS < 


oc Tau 


uooepyy 


u 


2iy.04io 




neon 1 1 Q c 
0530+ 135 


0530+135 




MWL> 34y/Mars 


OT Q 1 

27.34+8.12 


Q Tau 


05Scp99 


D 


Both 












25.78+5.64 


a Boo 


OlJuilOO 


D 


216.0982 


2.91 hr 


3C273 


3C273 


MWC 349 


Mars 


81.55+2.35 


a Boo 


OlJunOO 


D 


219.5418 


2.91 hr 


3C273 


3C273 


MWC 349 


Mars 


85.45+2.44 


a Boo 


OlJunOO 


D 


Both 












83.50+1.71 


a Tau 


09Dec00 


C 


106.6845 


3.19 hr 


3C111 


3C111 


W30H, Mars, 
Uranus 




12.55+2.16 


a Tau 


09Dec00 


C 


110.1282 


3.19 hr 


3C111 


3C111 


W30H, Mars, 
Uranus 




12.77+2.95 


a Tau 


09Dec00 


C 


Both 












12.40+1.91 


a Tau 


All 


All 


108.40 












13.97+1.46 


a Tau 


All 


All 


217.82 












25.78+5.64 


a Boo 


All 


All 


108.40 












20.09+0.69 


a Boo 


All 


All 


217.82 












83.50+1.71 



Table 2. ISOPHOT measurements of sky background near a Tau 



Filter 


A 


TDT 


Meas.# 


Background+Unc. 




(/im) 






(MJy sr-l) 


Cl-100 


100 


86401601 


2 


42.32+3.57 


Cl-100 


100 


86401602 


2 


44.53+2.68 


C2-160 


170 


86002102 


1 


43.06+4.39 


Cl-60'' 


60 


84300501 


5 


34.54+4.68 


Cl-100'' 


100 


84300501 


2 


36.69+4.14 


CI-6O'' 


60 


84300501 


5 


40.89+5.53 


Cl-60 


60 


combined 




37.19+3.57 


Cl-100 


100 


combined 




42.24+1.90 



'»P32 

''scaled P32 
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Table 3. ISOPHOT measurements of sky backgroTind near a Boo 



Filter 


A 


TDT 


Meas.# 


BackgrounditUnc. 




(fim) 






(MJy sr-l) 


Cl-50 


50 


27501511 


1 


18.16±1.42 


Cl-90 


90 


27502117 


1 


14.62±0.87 


Cl-105'' 


105 


27500602 


1 


12.54±0.54 


Cl-105 


105 


27500602 


1 


11.56±0.94 


C2-120 


120 


27503008 


1 


9.10±1.16 


C2-135 


135 


27503311 


1 


8.73±1.12 


C2-160 


170 


27503614 


1 


9.23±0.96 


C2-160'' 


170 


27503614 


1 


8.34±0.27 


C2-200 


200 


27502402 


1 


6.61±0.39 


Cl-105 


105 


combined 




12.30±0.47 


C2-160 


170 


combined 




8.41±0.26 



"■Central pixel only 



Table 4. Comparison of ISOPHOT sky backgrounds and fixed dark currents in W cm ^ fim 



Star 


ISOPHOT 


Nearest 


Fixed 


ISOPHOT 


DIRBE 




Band 


LWS det. 


Dark 


Sky 


zodi 


a Tau 


Cl-60 


SW2 


4.2E-18 


4.5E-19 


3.8E-19 


a Tau 


CI- 100 


LWl 


8.0E-19 


1.3E-19 


l.lE-19 


a Tau 


C2-160 


LW4 


9.1E-20 


4.5E-20 


5.9E-20 


a Boo 


Cl-50 


SW2 


4.2E-18 


2.0E-19 


>1.3E-19 


a Boo 


Cl-90 


SW5 


3.2E-18 


8.7E-20 




a Boo 


Cl-105 


LWl 


8.0E-19 


5.2E-20 


2.7E-20 


a Boo 


C2-120 


LW2 


7.7E-21 


1.2E-20 


1.2E-20 


a Boo 


C2-135 


LW3 


2.0E-20 


1.7E-20 


7.0E-21 


a Boo 


C2-160 


LW4 


9.1E-20 


1.5E-20 


1.8E-20 


a Boo 


C2-200 


LW5 


4.4E-19 


8.6E-21 
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Table 5. Radiative models: fluxes and brightness temperatures 



Frequency 


Quantity 


a Boo 


a Tau 


217.820 GHz (1.376 mm) 


Observed Flux (mjy) 
Predicted Flux (mJy) 
1 Predicted-Observed Flux | / c (observed) 


83.50±1.71 
33.43 
29.3 


25.78d=5.64 
32.38 
1.17 




Observed Brightness Temperature (K) 
Predicted Brightness Temperature (K) 


7040±140 
2820 


2200±480 
2760 


108.406 GHz (2.765 mm) 


Observed Flux (mJy) 
Predicted Flux (mJy) 
1 Predicted— Observed flux | /ct (observed) 


20.09±0.69 

7.94 
17.6 


13.97±1.46 

7.80 
4.23 




Observed Brightness Temperature (K) 
Predicted Brightness Temperature (K) 


6840±230 
2700 


4810±500 
2690 



Table 6. Radiative Models: atmospheric extents 



Quantity 


a Boo 


a. Tau 


Adopted angular diameter (mas) 


21.0±0.2 


20.88±0.10 


Adopted parallax (mas) 


88.85±0.74 


50.09±0.95 


Deduced stellar radius (cm) 


1.77 X 10^^ 


3.12 X 10i2 


1.64 (Um - 2.8 mm atmospheric thickness (cm) 


4.38 X 10^0 


4.22 X 10^0 


Percentage of stellar radius 


2.5% 


1.4% 



